Lattice solitons in quasicondensates 
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We analyze finite temperature effects in the generation of bright solitons in condensates in optical 
lattices. We show that even in the presence of strong phase fluctuations solitonic structures with 
well defined phase profile can be created. We propose a novel family of variational functions which 
describe well the properties of these solitons and account for the nonlinear effects in the band 
structure. We discuss also the mobility and collisions of these localized wave packets. 



Bose-Einstein Condensates (BEC) in optical lattices 
(OL) are unique candidates to explore phenomena that 
are often extremelly elusive in other areas of physics. The 
dynamics of the system is dominated by the interplay 
between the nonlinearity, that can be modified via Fes- 
hbach resonances [lj, and the periodicity which can be 
engineered through the intensity, geometry, polarization 
and phase of the two counterpropagating laser beams 
conforming the optical lattice. BEC in optical lattices 
exhibits phenomena known from solid state physics as 
demonstrated for instance in Q. Furthermore, its simi- 
larity to other cubic nonlinear periodic media has stimu- 
lated a renewed interest in solitonic models [H, 0, IS El Q • 

Solitons, which, strictly speaking, are exact solutions 
of integrable models corresponding to wave packets that 
propagate without change of their shapes and veloci- 
ties even in the presence of collisions, appear in many 
branches of physics. In one-dimensional homogeneous 
condensates with attractive interactions, bright solitons 
exist ||, and have been recently observed 0. The pres- 
ence of a periodic potential, like, e.g., an optical lattice, 
breaks the translational invariance, and the system be- 
comes most likely non-integrable, having less conserved 
quantities than degrees of freedom. Nonintegrable sys- 
tems, however, admit also localized solutions which are 
as well commonly termed solitons. These structures dif- 
fer from proper solitons either in their motion and/or 
their collisions. A well known example are optical soli- 
tons in periodic media (cf. [10]) whose interactions have 
been extensively studied using either the discrete non- 
linear Schrodinger (DNLS) equation [llj ] or weakly per- 
turbed integrable models |l2j . 

So far, to our knowledge, the generation of lattice soli- 
tons in condensates has been only discussed at zero tem- 
perature 0,0, IS 00) where there is an analogy between 
an array of optical waveguides in a Kerr medium and a 
BEC in a periodic potential. In this limit, one can gen- 
erate bright lattice solitons in repulsive condensates 
when the tunneling rate balances the nonlinear energy of 
the system. This compensation occurs if the soliton is 
placed at the edge of the first Brillouin zone where the 
effective mass becomes negative reaching the staggered 
configuration Q. 

Here we discuss the generation of solitons in repul- 
sive condensates at finite temperature. We show that 



even when the condensate is not phase coherent, solitonic 
structures with well defined phase profile can be created. 
A new insight into the nature of these lattice solitons 
is obtained by means of a novel variational ansatz that 
accounts for the effects of the nonlinearity. Finally, we 
address the issue of mobility and collisions. 

For BEC in 3D trapping geometries, fluctuations of 
density and phase are only important in a narrow temper- 
ature range near the BEC transition temperature T c |l5j . 
For pure ID systems, however, phase fluctuations are 
present at temperatures far below the degeneracy tem- 
perature Td = N hujjr /kn, while density fluctuations are 
still suppressed |ljj (ks denotes the Boltzmann con- 
stant, N the number of atoms, io x the axial trapping 
frequency). Phase fluctuations can be studied by solving 
the Bogoliubov-de Gennes (BdG) equations describing 
elementary excitations. Writing the quantum field oper- 
ator as ipix) — y/ no(x) exp(z0(x)) where no(x) denotes 
the density, the phase operator takes the form |l7|: 
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where a,j is the annihilation operator of the exci- 



and energy ej = 
<j, where Uj and Vj 



tation with quantum number j 
ftu xy /j(j + l)/2, and ff = Uj + < 
denote the excitation functions determined by the BdG 
equations. In ID and in the Thomas-Fermi (TF) regime, 
the functions ff have the form: 
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where Pj(x/ Rtf) are Legendre polynomials, Rtf — 
fin/mujx 2 ) 1 / 2 is the TF radius and \i the chemical poten- 
tial. The phase coherence length, = RTFTdhu> x / ^T, 
characterizes the maximal distance between two phase 
correlated points in the condensate. Phase fluctuations 
increase for large trap aspect ratios, u>t/u> x , and small N 

To study temperature effects in the generation of 
lattice solitons, we consider a 87 Rb condensate with 
N = 500 atoms in a magnetic trap with frequencies 
uj t = 715 x 2-7T Hz and lu x = 14 x 2tt Hz. For such pa- 
rameters, the system is effectively in a ID regime since 
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transverse excitations are supressed (fi « tuvt)- More- 
over, along the axial direction, the condensate is well in 
the TF regime (n » fuv x ). To create lattice solitons in 
repulsive condensates, the optical lattice depth belongs 
to the so-called weak potential regime. Thus, one cannot 
use the tight binding approximation to rewrite the 
condensate order parameter as a sum of wavefunctions 
localized in each well of the lattice. To study the dynam- 
ics of the system we use the ID Gross Pitaevskii equation 
(GPE): 



ihip(x, t) 
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with coupling constant g = 2Nha s L0t and scattering 
length a s =5.8 nm. The confining potential V(x, t) — 
mu> x {t) 2 x 2 /2 + Vo(t) sm 2 (irx/d) describes both the axial 
magnetic trap and the optical lattice. The later is charac- 
terized by its maximal depth Vb = lE r and by its spatial 
period d = A/2 (A = 795nm). Energy is expressed in 
units of the recoil energy E r — H 2 k 2 /2m, where k = w/d. 
Temperature is included at the level of the GPE by cal- 
culating first the density at T = in the presence of the 
magnetic trap only, i.e., V(x, t = Q) = mu> 2 x 2 /2 and then 
the phase operator (1). To this aim, we calculate the Bose 



ocuppation Nj 



1] 1 of the low-energy modes 



in the Bogoliubov approximation for a fixed T replacing 
the operators dj and a+ by random complex variables ctj 
and a* such that (|aj| 2 ) = Nj 



[l8| . One can also study 



the quasi ID case using a 3D GPE and replacing in (1) 
the energy of low excitations by Cj = hj x yj j(j + 3)/4 
[jfij and the Legendre by Jacobi polynomials (18j |. 

We summarize the procedure to generate lattice soli- 
tons at a given T. Once the phase fluctuations have 
been included, periodicity is introduced by growing adi- 
abatically the optical lattice (t< r). During this process, 
density fluctuations arise. We simultaneously turn off 
the magnetic trap and imprint a phase difference of tt 
between consecutive wells to place the condensate at the 
edge of the first Brillouin zone. This is achieved by us- 
ing an auxiliary optical lattice with double spatial period 
than the first one, acting for a time much shorter than 
the tunneling and the correlation time of the system Q . 
Hence, the condensate acquires the desired phase without 
modifying its density profile. 

Figure ^ shows density and phase of the condensate at 
T = 0.8T C (T c = Nhaj x /k B ln2N [H|) (a) in the pres- 
ence of the magnetic trap only (i = 0), (b) after growing 
adiabatically the optical lattice (t = r), and (c) 100ms 
after the phase imprint has been performed and the mag- 
netic trap has been switched off (t = t + 100). Figure 
1(d) shows a lattice soliton generated at T — with oth- 
erwise identical parameters. Despite the strong fluctua- 
tions of phase and density induced by finite temperature, 
the system (after phase imprint) evolves towards a stag- 
gered soliton configuration. It contains approximately 
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Figure 1: Density and phase (inset) profile of (a) magnetically 
trapped (u x = 14 x 2n Hz) 87 Rb ground state condensate with 
N — 500 at T = 0.8T C , (b) after the adiabatic growing of 
the optical lattice (Vo = 1E T and A = 795nm), (c) of lattice 
soliton 100ms after the imprinting of a phase difference of 
7r between consecutive wells and after the magnetic trap is 
switched off (d)of lattice soliton generated under the same 
conditions of (c) but at T — 0. 



35% of the initial atoms and remains so for times much 
larger than the tunneling time. This robustness can be 
understood by realizing that the size of the lattice soliton 
is smaller than the phase coherence length for all T <T C , 
The soliton size (which depends on N, being smaller as 
N increases ) is independent of the temperature but due 
to the random character of the fluctuations the position 
of the generated soliton is different for each realization. 

An analytical description of lattice solitons in repul- 
sive condensates is often performed through an effective 
theory for the soliton's envelope in the effective mass ap- 
proximation [2^] • The influence of the periodic potential 
is included there via an effective mass and coupling con- 
stant. Instead let us consider a very general ansatz for 
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the soliton wavefunction: 



ip(x, t) = e- iflt / h G(A, a, x, x )F(x) 



(4) 



where G(A, a, x, Xq) describes the envelope of a soliton 
of amplitude A, width a, and centered at xq. The ef- 
fect of the periodic potential is included, in the weak 
approximation limit, as a combination of linear Bloch 
functions F(x) — ^ k e, lkx fk{%), where the fk{x) have 
the periodicity of the lattice. At the edge of the first 
Brillouin zone, this combination can be approximated 
by only two harmonics, i.e., F(x) — cos(27rx/A) ||. 
Notice that such an ansatz does not take into account 
nonlinear effects. Moreover, it does not present a mini- 
mum in the energy functional of the GPE for any atom 
number. Inspection of Fig^c) and (d) shows that the 
density profile inside each well is shifted with respect to 
the minimum of the optical potential. A relative sim- 
ple function that reproduces such site dependent shift 
is F(x) — cos(27rx/A'(x)) with an effective wavelength 
A'(x) = A(l + a(A)|G(A,cr,x,x )| 2 ) where a ( N ) « 1 
is a variational parameter which depends only on N. In 
other words, the modification of the band structure due 
to the nonlinearity results in an effective change of the 
periodicity of the system. The ansatz (@J with the ef- 
fective wavelength X'(x) does not allow a fully analytical 
treatment for the energy functional. As a first approxi- 
mation we assume thus that the effect of the nonlinearity 
is to shift A simply by a constant A' = A + 6. With this 
new ansatz, imposing the normalization of the wavefunc- 
tion (i.e., conservation of N) and assuming a gaussian 
envelope G(A, a, x, xo) = Aexp((x — xq) 2 /2a 2 ) the en- 
ergy functional E[A, a, xo, S] reads: 
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where B = \A\ 2 ^a/A, k' = 2tt/X' and k± = k' ± k. 
We minimize (JjJJ for a soliton centered in one lattice site 
i.e., xo = 0, as a function of the shift 6 and the width a. 
The other free parameter A, the amplitude, is fixed by 
normalization, i.e., 1 = 2_B(1 + e ~ fe °" cos(2fe'xo)). 

For each TV and xo, minimization of E[A, a, xo, S] with 
respect to a and 5 presents now a clear minimum. The 
values of a and 6 obtained through minimization coincide 
with the values of the width and on-site shift obtained by 
solving numerically the GPE J2Jl. Also, minimization of 
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Figure 2: Collisions between two identical lattice solitons 
(N — 200) at an initial distance of 2xo = 246d and with an 
initial phase difference of (a) A(j> = Orad, (b) A(f> = 7r/2rad, 
(c) A(p = 0.3rad and (d) A<j> = 0.9rad. The initial transfer of 
momentum is O.lfc/t in all the cases. 



(JSJ correctly indicates that solitons containing few atoms 
(N < 100) are very extended (large a) and present prac- 
tically no shift (S < 0.001), as one should expect in the 
continuum limit. On the contrary for larger A, a de- 
creases and S grows. 

The Peierls-Nabarro (PN) barrier H H3| i.e., the en- 
ergy difference between a solitonic configuration cen- 
tered in one minimum of the lattice, E[A, a, xo = rid, 5], 
and the one centered in one maximum, E[A',a',xo — 
nd/2,6'], has been also calculated using (J5J. For any 
value of A, E[A, a, x = nd, 5] < E[A', a', x = nd/2, 5']. 
Thus, even the constant shift approximation reproduces 
the main features of lattice solitons. 

On-site lattice solitons are created initially at rest. 
Nevertheless, if they are well localized in momentum 
space and an instantaneous transfer of momentum (large 
enough to overcome the PN barrier) is given to the sys- 
tem, the solitons move in opposite direction to the given 
momentum due to their negative effective mass. No- 
tice that the variational ansatz F(x) — cos(2irx/\') with 
A' = A + S is meaningful only in the static case. Such 
a simplified solution cannot be used to study soliton dy- 
namics, since the soliton is always chirped with respect 
its center. The Euler-Lagrange equations for an ansatz 
whose periodicity depends also on the density are quite 
complex. Therefore, to study dynamical behaviour one 
has to rely on numerical simulations. 

Giving momentum to the system is generally accom- 
panied by radiation of atoms. The larger the given mo- 
mentum is the larger the losses. For small initial transfer 
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of momentum the system exhibits a nonlinear response, 
its velocity slows down, and eventually a complete halt 
of the soliton occurs. The given momentum p, for which 
the soliton experiences a nonlinear response lies on the 
range < p < h(k — k') where k' corresponds to the 
inverse of the effective wavelength. Thus the linear re- 
sponse is recovered for broad solitons since k ~ k' and 
for well localized solitons if the initial given momentum 
is large enough, i.e. p > h{k — k'). Recently nonlin- 
ear movement of discrete solitons has also been reported 
P| . The interaction between lattice solitons depends also 
strongly on the initial given momentum. We study nu- 
merically collisions between two identical solitons with 
N = 200 each, created initially at rest and separated 
by a distance 2xq = 246d For small values of the ini- 
tial momentum p ~ Q.lkh, the two initial solitons merge 
at the interaction point resulting in a single soliton (with 
N = 200) at rest. The excess of atoms is lost by radiation 
in a symmetric way (Fig. 2(a)). For intermediate values 
of the initial kick, p ~ 0.2kh, the two solitons merge in a 
moving soliton with smaller N than the initial ones. For 
larger kicks, p ~ 0.3kh, the solitons exhibit a quasi elas- 
tic collision where the two solitons pass each other. The 
average phase difference between the two solitons also af- 
fects the nature of their interactions. For p = O.lkH and 
identical initial phase distribution one obtains the fusion 
of the initial two solitons into a single one at rest (Fig. 
2(a)). This merging behaviour is always present for aver- 
age phase difference between the solitons below 0.19rad, 
although the final structure can either be a soliton or 
a breather and its position can also change. For phase 
differences of the order of tt/2 or larger, the two soli- 
tons always repel each other (Fig. 2(b)). Between these 
two extreme cases, the dynamics of the system is un- 
predictable. Two solitons either with approximately the 
same iV and same velocity or with different N (being the 
one with more atoms the slower) (Fig. 2(c)) are examples 
of possible situations after the collision. In other cases, 
the two initial solitons form a bound state that ends up 
into a single soliton (Fig. 2(d)). These results evidence 
the strongly dependence on the initial conditions in the 
dynamics of the system. 

Summarizing, finite temperature effects in the genera- 
tion of lattice solitons in quasicondensates have been ad- 
dressed. New insight into the nature of these structures 
is provided through a new family of variational functions 
in which the effect of the nonlinearity is shown as an ef- 
fective change in the lattice periodicity. We hope that 
such novel ansatz opens new possibilities in the study of 
lattice solitons. 
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